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We establish the most general form of the discrete elasticity of a 2D triangular lattice 
embedded in three dimensions, taking into account up to next-nearest neighbour 
interactions. Besides crystalline system, this is relevant to biological physics (e.g., 
red blood cell cytoskeleton) and soft matter (e.g., percolating gels, etc.). In order 
to correctly impose the rotational invariance of the bulk terms, it turns out to be 
necessary to take into account explicitly the elasticity associated with the vertices 
located at the edges of the lattice. We find that some terms that were suspected in 
the litterature to violate rotational symmetry are in fact admissible. 
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1. Introduction 

Since Born and von Karman, the construction of microscopic elastic models of 
crystals has formed an important part of solid-state physics, aimed to determine 
elastic constants and lattice dynamical properties (Born & von Karman 1912; Born 
1914). Imposing global translational and rotational invariance is a delicate matter, 
although the required relationship are well established (Born & Huang 1954). In 
fact, inconsistencies associated with rotational invariance in the earlier models were 
pointed out by Lax (1965) and Keating (1966a, 6). More recently, the importance 
of triangular lattices has emerged with the study of soft matter systems such as 
percolating gels (Feng & Sen 1984; Schwartz et al. 1985; Arbabi & Sahimi 1993), 
membrane with crystalline order (Kantor & Nelson 1987; Seung & Nelson 1988; 
Gompper & KroU 1997), biological membranes (Discher et al. 1997; Coughlin & 
Stamenovic 2003), red blood cell cytoskeleton (Saxton 1990; Lim et al. 2002; Mar- 
celli et al. 2005), and finite element modeling (Gusev 2004). 

In this paper, we consider a triangular lattice embedded in three dimensions and 
we determine the most general expression of the harmonic elastic energy up to next- 
nearest neighbour interactions. We formulate a systematic procedure to generate 
all the energy terms allowed by symmetry. This procedure could be applied to other 
type of lattices. We show that in order to correctly ensure rotational invariance, 
it is necessary to consider the boundary of the lattice, as first suggested by Lax 
(1965). Our central result, equation (3.2), gives all the allowed energy terms, with 
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no particular connection to mechanical elements such as springs, elastic wedges, 
etc. This contrasts with the above cited papers that postulate energy terms in a 
non exhaustive way. In particular, the most general energy contains a non-central 
nearest neighbour term (with coefficient K^), which was believed to violate rota- 
tional invariance in percolating gels (Feng & Sen 1984) on the grounds of Keating's 
argumentation (1966a, b). The situation is actually delicate since in a later paper, 
Keating (1968) restricted the generality of some assertions of its earlier papers. 

Our paper is organised as follows. In §2 we formulate the elastic energy of 
triangular lattices. In §2 a, wc write down all the scalar energy terms quadratic in 
the displacements of the lattice vertices. In §2 b, we retain only the terms satisfying 
the six-fold symmetry of an infinite triangular lattice. In §2 c, we consider also the 
edges: we study the symmetries relative to them and we retain the allowed edge 
terms. In §2d, we impose global translational and rotational invariance, taking 
explicitely into account the reduced symmetry of the boundary, and we retain the 
allowed terms. In §2e, we collect all the contributions and wc extract the resulting 
bulk terms. In §3, we analyse and discuss the results obtained in §2. In §3 a, we give 
a possible interpretation of our energy in terms of mechanical elements. In §3 b, we 
rewrite the energy into the generic expression (3.2), which constitutes our central 
result. In §3 c, we establish the connection between our model and its continuum 
version. Finally, in §4, we give a brief summary and we detail the origin of the 
rotational invariance problem. 

2. Discrete elasticity model 

(a) General definitions 

Let us consider a bounded triangular lattice. In its undistorted state, we assume 
that it lies in a plane V and that its vertices (i.e. its sites) are evenly separated by 

a distance £ (figure 1). We denote by r" € P the position of a generic vertex a. 
We also assume that the elastic properties of the lattice are symmetric with respect 
to v. 

In an arbitrary distorted state, r" becomes r" + it", where is a three- 
dimensional displacement. The most general quadratic expansion of the lattice 
elastic energy, F, is given by 

ir({u«}) = ^ti«.S«.u«-F^^tx«- (2.1) 

a a /3jta 

In this expression, we have separated the terms coupling identical vertices from 
those coupling different vertices. We have excluded the terms linear in it", because 
they must vanish for the undistorted state to be a minimum of the energy. Since 
F is scalar, the tensors 5" and A^f^ must satisfy (S")* = and (A"'')* = 
where the superscript t means transpose. Note that in general A""^ ^ A^". 

We fix the interaction range by assuming that all the A^^^ vanish except when 
a and /3 are nearest or next-nearest neighbours. 

(6) Bulk terms 

To begin with, we consider the bulk vertices, i.e. the vertices that are far away 
enough from the boundary (this point will be precised in §2 c) so that the triangular 
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symmetry may be assumed. In the bulk the elastic properties are assumed to be 
spatially homogeneous. Let n be a unitary vector normal to V. Let (ti, be two 
three-dimensional vectors forming a local orthonormal basis in V, allowed to vary 
from site to site. The tensors S" and A"^^ can be decomposed into the complete 
basis B = (ti ti, ti (g) t2, <2 (8) <i, t2 n n), where the symbol (8) represent 
the tensorial product. Elements linear in n, such as ti (g) n, have been excluded, 
since the imdistortcd lattice is symmetric with respect to V. 

To construct S" , wc take ti along any one of the six nearest-neighbour directions 
of the imdistorted lattice, and t2 = nxti. Because of the triangular symmetry, each 
term in the decomposition of S" must be even in ti and in <2 , and the coefRcients of 
ti (g) ti and t2 (g t2 must be equal (this condition is necessary for S°' to be invariant 
under an in-plane | rotation). Hence, 

= si Ja + sm O n, (2.2) 

where I2 = ti ti + t2 ® <2 is the identity tensor in the plane V. By virtue of the 
lattice homogeneity, the coefficients si and Si should not depend on a. 

We now construct A"^ . Let t"*^ be the unit vector in the direction going from 
to r^. We take ti = f^^ and t2 = n x ti. Without loss of generality the 
decomposition of A"''^ over B can be rewritten as 

j^a0 ^ ^ i^ + ^ri(E)n+^ ti(^ti + ^ (ti(gt2+<2(S)ti) + ^(*i«)*2-*2<8*i)- (2.3) 

By virtue of the lattice homogeneity, the coefficients {a^, a^, . . .} take only two sets 
of values: either {oi, ai, . . .} if a and /3 are nearest neighbours, or {02, 0,2-, . . .} if a 
and P arc next-nearest neighbours. The in-plane mirror symmetry with respect to 
requires that A"-^ must be even in ±2, which leads to 

5i = a* = 0. (2.4) 



(i) Collecting bulk terms 

We now anticipate that we shall clearly distinguish between the S" and A"^ 
belonging to the bulk and those belonging to the edges. Denoting by Fbuik the part 
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Figure 2. The seven types of neighbour clusters as indicated by the large dots (the lines 

are only guides for the eyes). 

of the energy F which involves only the bulk terms, we arrive, from equation (2.1), 
(2.2), (2.3) and (2.4), at 

(bulk) 

a 

(biilk) 2 
<al3> i=l 

where = • n denotes the normal component of u" and = — u"n the 
projection of u" on P. The symbol J2<a/3> indicates the summation over both 
nearest and next-nearest neighbour pairs of vertices. As already mentioned, we 
assign to the former the suffix i = 1 and to the latter the suffix i = 2 (see figure 2 
where these numbers are indicated as I and II). The function 5"^ takes the value 
1 if the pair {a, /?} is of type i and zero otherwise. We shall hereafter call such 
functions indicatrix functions. 

(c) Edge terms 

With arbitrary coefficients si,si,ai,ai and Si {i = 1,2), the energy form (2.5) 
is not invariant under a global translation or rotation of the lattice. In order to 
correctly ensure rotational invariancc, it is necessary to take into account the edges 
of the lattice. This seems counterintuitive, because one expects the contribution of 
the edges to be subdominant with respect to the expensive contribution of the bulk. 
However, as it will more precisely be shown in §4, this is not the case, because in a 
rotation the displacements of the edges are proportional to the size of the system. 

(i) Neighbour clusters and vertices types 

Since our interest lies in the bulk terms, we choose the simple shape of a regular 
hexagon of edge length N£ (figure 1) as our system. Consistently with the energy 
expansion, the specificity of the edges must be taken into account up to the next- 
nearest neighbour distance. For this purpose, we group the vertices into neighbour 



:ul-ui + ai < + di {ul ■ f^") (< • t"^) , (2.5) 
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Figure 3. The 5 types of vertices (left) and the 8 types of vertices doublets (right) together 
with their local basis ti, t2 (the former being in thick arrows). Only the basis relative to 
the edges are indicated. For the doublet (a,/3) of types 3,4,7 and 8, the begininng (resp. 
end) of the thick arrow corresponds to the first (resp. second) component a (resp. 0) of 
the doublet. 

clusters in which aU the intcr-distanccs arc at most equal to the next-nearest neigh- 
bour distance {t^/2>). There are seven types of such clusters, which we label from I 
to VII (see figure 2). A vertex is said to participate in a neighbour cluster of type 
Y if it can be grouped with its neighbour vertices into a neighbour cluster of type 
Y. Note that vertices near the edges participate in a smaller number of neighbour 
clusters than those in the bulk, because of their lacking neighbours. 

To each vertex a we now assign a type depending on which neighbour clusters it 
participates in: two vertices are assigned the same type if, and only if, their positions 
within the set of neighbour cluster in which they participate are indistinguishable. 
Close inspection shows that the edge vertices are classified into five types, labeled 
from 1 to 5, as shown in figure 3. The vertices of type 1 are the bulk vertices. 

We apply a similar procedure to assign a type to the oriented nearest and next- 
nearest neighbour vertex pairs, which we call vertex doublets. We distinguish the 
vertex doublets (a, j3) from (/3, a) since generally A"^!^ ^ A^^". Two vertex doublets 
are assigned the same type if, and only if, their positions and orientations within 
the set of neighbour clusters in which they participate arc indistinguishable. Wc 
thus determine eight types of vertex doublets, labeled from 1 to 8 (figure 3). The 
vertex doublets of type 1 are the bulk nearest-neighbour vertex doublets and those 
of type 2 are the bulk next-nearest neighbour vertex doublets. We also define types 
for the vertex pairs in the following way. When the two doublets associated with a 
given pair are of the same type, the pair is assigned that type (this holds for the 
types 1, 2, 5, 6) while the pairs corresponding to doublet types 3 and 7 are assigned 
the type 3 and those corresponding to doublet types 4 and 8 are assigned the type 
4. 



(ii) Symmetry allowed edge terms 

In order to take into account the local symmetry up to the next-nearest neigh- 
bour distance, wc consider the symmetries of the set of neighbour clusters in which 
the vertices or vertex doublets participate. Let us first examine the coefficients of 
the tensors for vertices of type 2 to 5 (edge vertices). Since each one possesses an 
axis of symmetry, we choose ti along this axis and -L ti as indicated in figure 3. 
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Then, as S"' must be even in t2, we arrive at 

S°' = Sil2+S,n(g)n + s^T°' (g>T", (i = 2,...,5) (2.6) 

where r" may either ti or t2- For definiteness, we choose r" = t2 for vertices of 
type 2 or 4 and r" = ti for vertices of type 3 or 5. 

Wo now construct the tensor A"^ for the vertex doublets of type 3 to 8 (edge 
types). The most general expression for A"^ is given by an expression similar to 
(2.3): 

+ ^(tiOt2 + t2Oti) + -|^(ti0t2-t2Oti). (2.7) 

We choose ti = t"^ and t2 as indicated in figure 3. The vertex doublets of types 4 
and 8 have a mirror symmetry with respect to ti, hence A°'^ must be even in t2, 

which implies = 0. Other relationships can be obtained by comparing 

the coefficients of A'^f' and those of A^". First, note that {A°'^y = A''" implies 

the vertex doublets of type 5 and 6 do not change type when exchanging a and /?, 
their coefficients must satisfy a"^ = a^", a"^ = a^", a"^ = a^"-, a"^ = 0.^°" and 
a*«^ = a*^". Therefore a"^ = for types 5 and 6. Finally, we find 

A"'5 = |j2 + |n®n + |r'3®r'3 + ^e"^^ + |a"^ (z = 3,...,8) (2.8) 

where e"'' = (g) t2 — 12 <Xi *i , cr"^ = ti (g) t2 + *2 ^ *i and = = 0,4 = dg, = 0.5 = 
ae = 0. In the following, we define e as the tensor that performs a counterclockwise 
^ rotation in the plane V; depending on the relative orientations of ti and t2, e"^ 
is either equal to e or to —e. Note that in general, we do find A"^ ^ A^", contrary 
to what is the case in the bulk. 



(iii) Collecting edge terms 

The total elastic energy can now be written as = i^buik + ^cdgo, with -Fbuik 
given by (2.5) and Fedge obtained by collecting the edge terms defined above. We 
regroup the terms u" ■ A"^ ■ and ■ A^"- ■ u°' in order to sum on the pairs and 
not on the doublets. We arrive then at 



a i=2 



<al3> i=3 



<a/3> 



(5f ??"^ci3<-cr"'3-Mi+ E 5fa*ul-e^ 



=3,5,6 



(2.9) 



In this expression, 5f € {0, 1} and 5f e {0, 1} are the indicatrix functions of the 
vertices of type i and of the vertex pairs of type i, respectively (these types are 
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defined in (i)) , and the symbol indicates the summation over ah pairs 

{a, /?} of vertices. The third line in (2.9) requires more explanations. We define a 
closed contour by grouping the vertex pairs of types 3 (sec figure 3), another closed 
contour by grouping the vertex pairs of type 5, and three separate closed contours by 
grouping the vertex pairs of type 6. The symbol £ denotes the sums over all these 
contours followed counterclockwise. The symbol equal to 1 or — 1 according to 

whether the direction from a to /3 goes outwards or inwards, respectively. Finally, in 
order to simplify the notation, we have renamed |a4 + ^ag ^ 04, ^0.4+ ^as ^ 0,4, 



\a4 + lag 



+ ^ag 

0,4, ^flg + ia* — > and \a%, — ^a* — > at. 



2U.3 2"7 



^3. 



{d) Rotational and translational invariance 

Our task is now to ensure that the elastic energy remains unchanged if the 
lattice is translated or rotated as a whole. Although the general relationships are 
well-established (Born & Huang 1954), it is not an easy task to implement them in 
such a way as to put the elastic energy in a satisfactory form. Our strategy is to 
rearrange F in order to put forward a number of terms that are obviously invariant, 
then to apply deliberately chosen translations or rotations in order to show that 
the remaining terms should vanish. 

Wc first separate in F (= -Fbuik + -Pcdgo) the contribution depending on the 
components from the contribution F^-^^ depending on the components . We 
thus have F = F^") + F^-^) with 



a i=l <a/3> i=l 



a 1=1 
6 



<a/3> i=l 



(2.10) 



<a/3> 



=3,5,6 



■a/3 



(2.11) 



(i) Translation invariance 

Putting forward obviously translationally invariant terms in the form of differ- 
ences, F^"^ can be rewritten as 



^^"^ = -J E E'5r"«^«-<)'+EE'^"^»«)'' (2-12) 



<a/3> i=l 



where the Ci are linear combinations of the coefficients Sj and cij (the coefficients 
in (2.12) still form an independent set). In the appendix, §a, we show that the 
invariance of F under any translation along n requires: 



Vz, d = 0. 



(2.13) 
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Similarly, F^-^) can be written as obviously translationally invariant terms plus 
remainders: 



<Q/3> i=l 



- E (1 - Cv,6 
<a/37>iv 

+ E E^l'^n^-o-^-K-o 

<a/37>in i=3,5 

+ EE'^"[^*(«")'+^^«-^")']- 

a i=2 



(2.14) 



Whereas the first line was clearly obtained by putting forward squares and discard- 
ing the remainders in the last line, the other lines deserve detailed explanations. 
Here and afterwards, the symbol X]<a/37>Y denotes the sum over all the triangular 
neighbour clusters of type Y, with a, j3 and 7 taken in counterclockwise order. For 
type IV, the vertex (3 is supposed to be the central one, as in figure 2. The symbol 
Sy^P" is the indicatrix function of the neighbour clusters {a, /?, 7} of type Y lying 
within the closed contour made by the vertex pairs of type i (i = 3,5 or 6). More 
precisely for i = 6 it means to lie within at least one of the three closed contours 
made by the vertex pairs of type 6. In the following, wc shall often encounter the 
translationally and rotationally invariant expression (as in the second line): 



q(0) 
^a/37 



(< - O X t^^" - (w^ - mX) X +0{u'^), (2.15) 



where da/Sj and 6* 



(0) 



? or ^ are the angles made by the vertices triplet {a, (3, 7} 



with f3 being the apex, after and before an elastic distortion of the lattice, respec- 
tively. In the appendix, § b (i), wc show how the terms involving er"^ in (2.11) were 
rewritten with the help of (2.15) and we give the expressions of the coefficients 
bi and h that were modified by this procedure. Finally, the lines involving e were 
obtained as detailed in the appendix, §b (ii). The coefficients Di and Di are linear 
combinations of the coefficients Sj, s *, and 03 and the new coefficients 

in (2.14) still form an independent set. In the appendix, § b (iii), we show that the 
invariance of F under in-plane translations requires: 



Vi, A = and A = 0. 



(2.16) 



Consequently, the total energy can be written in the following form, obviously 
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invariant under any translation of the system: 



<a/3> i=l 

- ^ (1 - <76)4 ^3 [(«^ - Ul) X t'^" - (nl - <) X t/^^] ' 

<a/97>iv 

+ E E<i^n<-o-e-(<-o 

<a/37>m i=3,5 

+ E <6"^s«-o-e-(<-o 

= +i^(-L)^ (2.17) 



where F^"^ is the first term and i^(^) is the sum of all the other terms. 



(ii) Rotation invariance 

Putting forward rotationally invariant terms, F^"'\ as obtained from (2.17), can 
be rewritten as 



F^") = -fc 1 V 6»2 

<a/37(5>vi 



E ^"/St* 

<a/37(5>vn 



E E 

<a/3> 1=2,4,5,6 



+ 0(^3), 



(2.18) 



The symbol X^<;q/375>y denotes the sum over all the neighbour clusters of type 
Y=VI or VII with the vertices placed as shown in figure 2. The angle BaSi-^s denotes 
the wedge angle between the planes defined by the vertices (a, /3, 7) and the vertices 
(a, /3, S). In the appendix, § c (i), we show how the wedge angles are related to the 
normal displacements The coefficients Ei arc linear combinations of fij, fii and 
as, and the new coefficients still form an independent set. We show in the appendix, 
§ c (ii), that the invariance of F with respect to out-of-plane rotations requires: 



Vi, Ei = 0. 



(2.19) 



Finally, we put forward in F^-"-) seven rotationally invariant terms, obtained by 
considering various energy terms associated with deviations of lengths or angles 
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with respect to their equilibrium values: 



<al3> i=l <a/37>in 



c.p. 



-'a/37 



c.p. 



<a0j>ui 

<a/37>iv 



c.p. 



2 ^ 



<a/37>iv 



7r\ 



IV,6 
2 



277 



'a/37 



«7 g; 

<a/37>iv 

+ E E<lGi«-«l)-^-K-0' 

<a/37>in i=3,5 

+ E <e « -O-e- {ni - + G ^ 5f{nl - n^f . (2.20) 



<a/37>v 



<a/3> 



The abbreviation c. p. stands for 'circular permutations' of {a, /3, 7} and ia0 is the 

distance between the vertices a and /?. The detail of how (2.20) and (2.17) arc 
related is given in the appendix, § d (i). The coefficient G is a linear combination 
of 63 and the other 6i's and the coefficients G* are linear combinations of the 
corresponding b* and the &j's {j 7^ 3). The new coefficients still form an independent 
set. In the appendix, §d(ii), we show that the invariance of F with respect to in- 
plane rotations implies 



G = and V«, G* 



0. 



(2.21) 



(e) Final expression of the bulk terms 

From the preceeding section, it follows that the most general elastic energy F of 
an hexagonal triangular lattice, invariant with respect to translation and rotations, 
is given by the sum of (2.18) and (2.20), with all the £^i's, GJ"s and G equal to zero. 

Putting forward the bulk terms, and expressing them in terms of the displace- 
ments instead of angles and lengths (using (2.15) and the the formulas of the ap- 
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Figure 4. Interpretation of the elastic terms in (2.22) in terms of linear and angular springs: 
The top drawing is in perspective while the others are in projection. The arrows are normal 
to the facets. The springs with elastic modulus k (resp. k') have an equilibrium length 
£ (resp. £y/3). The angular springs with moduli k„i and fc„2 (elastic wedges) have zero 
equilibrium angles; the angular springs with moduli fc^ or k'^ have equilibrium angles i/;. 

pendix, §c(i) and §d(i)), we finally obtain 



<Q/37(5>vi 
<af3j5>vii 



u?, - ul 



<al3> 



1 



Q/3 

2 

£2 



2^ /2 "1 

<a/37> 



^2 "111,3 
1 



<Q/3> 



(mJ - u^) X - (tt-J - U 



1 2 



+ it' V — A"'''^ 

+ 2^2^ 3^2 '^V,6 



<Q/37>v 



(w^ - tt^) X - (m'^ - mX) X « 



a7 



1 2 



C.p. 



C.p. 



i% ^ <5fvl ^ - u^^) X (ul ul) X 



(2.22) 



</3c[7>iv 

where k — ki and fc' = k2 and the ellipsis represent terms involving only edge vertex 
pairs. 



3. Analysis 

(a) Comments on the energy form (2.22) 

As shown in figure 4, each term of (2.22) can be given a simple interpretation in 
terms of linear or angular springs. The terms with coefficients fc„i, corresponding 
to elastic wedges, tend to keep the lattice flat and yield a form of curvature energy. 
The terms with coefficients k and k' represent elastic bonds acting between nearest 
and next-nearest neighbours, respectively. The terms with moduli kj^ , fc'^r and fci^r 
are angular springs acting at the bond joints. 
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The decomposition (2.22) is in many respects non-unique. For instance, the term 
with coefficient kn2 could be omitted if we are only interested in the bulk terms, for 
it has the same decomposition on the bulk vertex pairs as the term with coefficient 
kni; this can be easily checked by comparing the coeSicients of 5"'^ and Jj'^ in the 
equations of the appendix, § c (i). Note, however, that the expressions of these two 
energies in terms of clusters of four vertices differ, as we can see in the two first 
terms of (2.22). As another example, we could have included in (2.22) energy terms 
associated with the | angles in the bulk neighbour clusters of type IV; this would 
bring nothing new, since they can be written as a linear combination of other terms 
that are already present: 



-'a/37 



</3a7>iv 



Nonetheless, although the form of (2.22) depends on the choices we made during 
the various transformations on the energy, our procedure guarantees that it does 
correspond to the most general elastic energy of a triangular lattice up to next- 
nearest neighbour interactions. 



(6) Compact expression of the elastic energy 

With the help of the formulas given in the appendix, § c (i) and § d (i) , the energy 
(2.22) can be rewritten in a more compact form, in terms of the normal, transversal 
and longitudinal components of the relative displacements between nearest and 
next-nearest neighbours : 

^ = ^e{^"«-<)'+^*[( 

+ edge terms, 

where t"^ = n x t"^, and the symbols J2<ai3>i ^^'^ S<a/3>n indicate the sum 
over all the neighbour clusters of type I and II, i.e. over the pairs of nearest and 
next-nearest vcirtices, respectively. The normal coefficient, which corresponds to 
curvature energy, is given by Kn = ^~^(4A;„i + 12A;„2)- The transverse coefficients 
Kt and are given by Kt = ^-2(|A;^/3 -|- 6^2^/3) and Kj. = i~'^{^k'^/3 - §^2^/3), 
and the longitudinal coefficients Kf and K'^, are given by = fr'^{k + ffc^/a) 
and K'^ = £"^{k' + \k'^j^ + \k2Tt/z)- These new coefficients are independent linear 
combinations of the former. Moreover, the five functions of (proportional to 

the coefficients Kn, Kt, Kn, K[ and K'^), appearing in the form (3.2), can be shown 
to be linearly independent in their functional space. Therefore, the seven (or more) 
rigidities which one could associate with the mechanical elements drawn in figure 4 
reduce to only five independent elastic constants. 



:<-<).tf 



\Ki [«-<)-r^'' 



{ui - u^) • tf] % K[ \{ui - u1_) ■ r^' ' 



(3.2) 
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A number of comments follow. The angular coefficients of (2.22), i.e. fc^^g 
and fc27r/3, are present both in the tranverse and the longitudinal coefficients: this 
is because in a triangular lattice one cannot change the angles between bonds with- 
out changing the distances between the vertices. Conversely, the coefficients k and 
k' appear only in the longitudinal coefficients. Also, as explained above, the two 
coefficients kni and fc„2 collapse into a unique normal coefficient Kn. Most impor- 
tantly, note that the non-central term with coefficients iff, which was understood 
from Keating (1966) to violate rotational invariance, is actually allowed (even if we 
restrain the model to nearest neighbour interaction only). We shall comment on 
this in §4. 

Finally, although (3.2) is expressed as a sum of squared terms, the stability issue 
is not trivial, since these terms arc not independent. For instance, the stability of the 
undistorted lattice does not require all the five coefficients of (3.2) to be positive. For 
example, if all the coefficients in (2.22) are positive, which is obviously a sufficient 
condition for stability, and if fca^L > ^k'^,^, we may have K'^ < 0. Further comments 
on the statibility issue will be given in g4. 



(c) Link with the continuous theory 

In a continuous description, the displacement of the point r on the plane V 
is described by a three-dimensional continuous function u{r). In order to link the 
discrete description to the continuous one, we assume that the displacements u" 
vary very slowly at the scale of the lattice spacing £ and we identify u°' and u{r°'). 
Hence, we transform (3.2) with the help of the Taylor expansion: 



u 







Vti(r")-hC'(e2), (3.3) 

where — r°' = ^^t"^ , with ^ = £ for nearest neighbour vertices and ^ = £^/^ for 
next-nearest neighbour vertices, and 2f ~^/%/3 / d^r. We take an orthonormal 

coordinate system (.x,y, z), with the z axis parallel to n and the x axis along one 
of the three nearest-neighbour directions of the undistorted lattice. Then, for the 
in-plane part F^-"-' of (3.2), we obtain F^-^) = /d^r/^^', with 



+ Kt 



3 {u%f + 3 (,4)' + {u%f + {uy^f + 2 u% + 2 u% 

«x) ' + («f,) ' + 3 («fj ' + 3 («f,) ^-2ul uy^ - 2 ul , (3.4) 



in which Ki = \\/2,(K(, + 3K'^) and Kt = |\/3(A't + SK'^) and w*^ is the derivative 
of the zth component of u with respect to the jth component of r (note that 
i G {x,y,z} whereas j G {x,y}). Using integrations by parts, we can rewrite the 
energy density in the standard, isotropic, form = ^XeuBjj + fiSij Cy, with 
Bij = |(wj + u^i). The Lame coefficients are given by 

A = ^[Ke + 3K'f-5iKt + 3Ki)], (3.5) 

M = ^[Ke + 3K'i + 3{Kt + 3Kl)]. (3.6) 

In terms of the spring constants of (2.22), this gives A = |\/3^~^(gfc-|- ^fc' — A;^/3 — 
^w/3~^27r/3) aiid/U = |\/3-f~^(gfc-|-ifc'-|-fc^/3+fc^^3+fc2,r/3)- For a System composcd 
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only of springs between nearest neighbours, we recover the classical relation A = 
H = \\/^k/P (Kantor et al. 1987). As expected, the area compressibility modulus 
Ka = A + /i depends only on the spring constants and not on the angular constants: 
A + At = \V2,{Ki + iK[ -Kt- 2,K[) = ^V^i~^{k + 3k'). 

For the normal part F^"^ of the energy, it is necessary to add the second order 
terms t"^ ■ Wm • t"^ in the Taylor expansion because the first order vanishes. 
With f = / d^r /^"\ we obtain the isotropic energy density 

f^"^ = l^K. + n%y)\ (3.7) 

with the bending constant n = Kn = ^{kni + 3fc„2). One recognises the 

bending energy associated with the curvature H = ^{u^^x + "fyy) + 0{u^) of the 
surface z = u^{r). 



4. Summciry and discussion 

In this paper, we have determined the most general form for the elastic energy of 

a triangular lattice when both first and second neighbour couplings are taken into 
account. We have found that, among the six coefficients originating from the de- 
composition of the energy in terms of the normal, transversal and longitudinal com- 
ponents of the relative displacements between nearest and next-nearest neighbours, 
only five are independent (equation (3.2)). Consequently, only five combinations 
of the elastic parameters given for instance in figure 4 are measurable. We have 
also found that only three coefficients survive in the large scale limit, the Lame 
coefficients and the bending modulus; therefore our discrete description introduces 
two new elastic coefficients. From the following well-known stability criteria: k > 0, 
/i > and A + /i > 0, it is possible to deduce some necessary stability conditions: 
Kn > 0, Ke + 3K'f > and Kt + SK'^ > -^K^ - K[. A complete stabifity analysis 
can be performed by Fourier transform. However, the additional stability conditions 
on the coefficients being rather intricate, we shall not describe them here. 

One of the main issue in the paper was the necessity to take into account the 
boundary elastic energy in order to correctly impose rotational invariance. This 
is well exemplified by the following energy term, which represents the sum of the 
energies of the angular springs between nearest bonds: 

E(».-j)' = i:l(<-<r 

+ Ei(<-<) +E (4-1) 

where i runs over all the pairs of adjacent bonds, making an angle 9i, and a—/3 runs 
over the interior bonds when (B) is indicated, or over the boundary bonds when 
(dB) is indicated. In the last sum, the bonds (a^P) must be taken counterclock- 
wise. Indeed, the left-hand side is manifestly rotationally invariant and, under an 
infinitesimal in-plane rotation u']_ = un x r", the last term of the right-hand side 
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Type 


Vertex 


Pairs of vertices 


1 


Vi =3N'^-m + 7 


Pi = gN'^ - 157V + 6 


2 


V2 = 6N- 12 


P2 = 9N'^ - 157V + 6 


3 


V3 = 6 


P3 = 12A^ - 12 


4 


^4 = 6iV - 6 


P4 = 6 


5 


^5 = 6 




6 




Pe = 12iV - 6 



Figure 5. Number of vertices, Vi, and of pair of vertices, Pi, as a function of their type 
(see figure 3) for a regular hexagon of edge length Ni. 

(r.-h. s.), which is a boundary term, can be shown to yield an extensive contribu- 
tion —2uj^ A, where A is the area of the lattice, just as the first term of the r.-h. s, 
which is a bulli term. It follows that: i) this first term is not rotationally invariant 
by itself, ii) that boundary terms are therefore necessary to correctly ensure rota- 
tional invariancc, and that, as announced, iii) the term with coefficient Ki is indeed 
allowed, even if we restrain the model to nearest neighbour interaction only. 

In conclusion, the method developed in this paper — which consists in defining 
neighbour clusters based on the interaction range, in considering explicitcly the 
edges of the lattice, in considering the local symmetries defined in terms of cluster 
environment to distinguish the various edge terms, etc. — provides a safe basis to 
further implement rotational invariancc. This method could be applied to a large 
variety of lattices and to higher dimensions. 

Appendix A. Translational and rotational invariances 

This appendix contains the detail of the calculations made to impose translational 
and rotational invariances to the elastic energy. 

(a) Translation parallel to n 

Here we demonstrate (2.13), on the basis of (2.12). Substituting u" = in 
equation (2.12) yields p(")/^^ = CiVi, where the coefficients Vi are defined 
in the table of figure 5. Hence, with ^ = 1, 

p{n) ^ 2.N'^Ci + 2.N (-9Ci + 6C2 + 6C4) + 7Ci - I2C2 + 6C3 - 6C4 + 6C5 (A 1) 

Setting P(") = for all N implies = Ci = C2 + C4 = C3 + C4 + C5. Next, 
imposing the same translation but from a distorted state, defined by it" = except 
for one vertex ao of type 4 (resp. 5) for which it"" = ^o"-, it is straightforward to 
show that the invariancc of F requires C4 = (resp. C5 = 0). Consequently, all the 
coefficients Ci must vanish. 

(6) Translation parallel to V 

Here we display in (i) and (ii) some details of the calculations made to go from 
(2.11) to (2.14) and we demonstrate in (iii) the equation (2.16) on the basis of 
equation (2.14). 
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(i) Treatment of the terms involving cr"^ . 

In the process to obtain equation (2.14), the terms involving cr"^ in (2.11) have 
been eliminated using the following identity: 

E (i - <6) 4 [i< - o X t'" - i< - o X 



<af5j>iv 



E(i-<6) 



<a/37>iv 



3\/3 



« - <)• tf^^+"— ]+{ul- <)• ( t^"+^— 



• 25"^' - 6"^ 



■'5 "^6 
2 



<a/3> 

[Sf - 2Sf + Isf^ - u^) . t"^]'^ 



(A 2) 



The coefficients bi and bi in equation (2.14) are then given by 6i = ai, 6i = ai, 62 = 
0,2, 62 = 02, &3 = a3-2\/3 53, 63 = a3 + 2\/3i3, 64 = 04, 64 = 04, 65 = as-^VsEa, 
ba = ns + 4a/3S3, fte = ae + \/353, fte = ae - |\/353, 63 = ag + S3, 65 = 05 + 2S3, 
6g = Og - as. 



(ii) Discrete Stokes-like formula: treatment of the terms involving e 

In the process of the transformation from (2.11) to (2.14), the terms involving 
e have been rewritten using the following general Stokes-like theorem transforming 
a sum of edge terms into a sum of bulk terms: 

£ dful-e-f'.u'i = J2 <,"(«i-0-e-(txl-0, (A3) 

<a/3> <a/37>Y 

where Y=III for i = 3, 5 and Y=V for i = 6, and where the vertices a, f3 and 7 are 
taken in counterclockwise order in the right-hand side of the equation. 



(iii) Invariance requirements 

Here and afterwards, t will denote any of the unitary vectors joining two nearest- 
neighbours vertices of the imdistortcd lattice. First we substitute u" = in equa- 
tion (2.14), which yields a third-order polynomial in N. Setting F = for all N 
implies = Di = 2D2 + D2 + 2D4 + 04 = 2D3 + D3 + 2D 4 + D4 + 2D5 + D5 , which 
leaves six independent coefficients among the Di and Di. Then, by imposing the 
same translation but on properly chosen distorted states, it is easy to show that 
all the coefficients Di and Di must be set equal to zero to satisfy the translational 
invariance. 
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(c) Rotation with the rotation axis parallel to V 

In this paragraph, we (i) display some details of the calculations made to go from 

(2.17) to (2.18) and (ii) demonstrate the equation (2.19) on the basis of equation 

(2.18) . 

(i) Decomposition of the terms involving wedge angles 
Expanding the first and second terms in (2.18) yields 

<a/37(5>vi <a/37<5>vi 

= ^ E {^^'^^ - + '^^f + ^'i^ + '^^f - -^e") « - <f+ 0{u% (A 4) 

<a/3> 

and 

E ^7*= E -^{2,U-n-<-ul-uif + 0{u^) 
<a/37i5>vii <a/37i5>vii 

= I E (S'^i" - 2'52" + + 3<5r^ - 5f) « - + 0{u^). (A 5) 

<a/3> 

The coefficients kni are then given by A;„i = |^^(ai — 203) and A;„2 = —^-^^(251 — 
Sas). 

(ii) Invariance requirements 

In an infinitesimal rotation with the rotation axis parallel to one of the directions 

of the undistorted lattice, wc have u" = 9t x r", which is parallel to n. Hence 
u°' — u^ = —9£txt°'^ for nearest neighbours, and u^ — u^ = —Ot^txf^^ for next- 
nearest neighbours. Substituting these expressions into (2.18) yields F^"^ / {^O^i'^) = 
3E2P2 + E4P4 + Er^Pr^ + SEePe, where the coefficients Pi are given in table 5. Setting 
F^"'> = for all N implies E2 = 0, E5 = -6Ee and E4 = SEq, which leaves only 
one independent coefficient. Next, imposing the same rotation but on a properly 
chosen distorted state, it can easily been shown that the last Ei must be set equal 
to zero to satisfy the rotational invariance. 

(d) Rotation within the plane V 

In this paragraph, we (i) display some details of the calculations made to go from 
(2.17) to (2.20) and (ii) demonstrate the equation (2.21) on the basis of equation 
(2.20). 

(i) Expansion of the various terms in (2.20) 

The first term in (2.20) is expanded as follows: 

6 6 , 

+ 0{u^). (A 6) 

<a/3> j=l <a/3> i=l 
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Setting A:27r/3 = 2&2 + 466 — 5^4 + 5^5)) the relations between the coefRcients 
ki and the coefScients in (2.17) are found to be ki = (?{—bi + &k2-K/i), ^2 = 

-£2(62 + 2^2^/3), fcs = f{-h + 3^2^/3), ki = e{-hi + 3fc2V3), ^5 = -£^&5, 

fee = -(."^{hQ + ^2^/3). 

The terms with coefficients /c^ , fc'x and k'i in (2.20) can be expanded using the 
identity (2.15) and the following identity: 



^ A"^T I [(< - u^^) X t''" - (m^ - O X t'^T] % c. p. I 



_ 3 
~ 2 



<a/3> 



<a/37>Y 



where either Y = III, A«^t = 5"^^J, c"^ = 5"^ + ^df,orY = Y, A«^t = Jyfg^, 

(-a/S — O2 -|- ; 01 1 — 111, LA — Ojjj 5 Ojjj 3, Ca0 — 2^3 T O4 • J-t 

follows that the relationships between the coefficients are = ("^{—^bi + 862 — 



1666 



b6 - 



ue -r I64 - §65), fc;/3 e^i-bi + 265 - 862 + 1266), k'^^^ = £2(462 
The terms with the coefficient k^^ in (2.20) can be expanded using the identity 
(2.15) and the following identity: 



1 2 



E A'^'^T [(< - ul) X t''" - (u^ - ul) X t'^T 

<a/37>iv 

= 3 E {d""^ (ul - 1x^)2 + d'"'^ [(u^[ - w^) • t"^ 



<a0> 



<a/37>v 



(A 8) 



where A'^/'t = Jjf^^, = 25^'^ - ^Sf + Sf + 6f - \5f, d'"/^ = -2Sf + 

2 , 1 <raf3 

'S"2 ~ °3 ~ "4 3"6 • 

The terms with coefficient fc'a^ in (2.20) can be expanded using the identities 
(2.15) and (A 2). It follows that fc^^/3 ^2^262 - 466 - ^^3). 

Finally, the term with coefficient fciL in (2.20) can be expanded using 



E (l-<"a) 



</3q7>iv 



+ 



(<-Oxt^"--^(<-nl)xt'^^ 
«-Oxt^"-i=«-<)xr/5 



1 2 



E <j e"'^ « - <)2 + e'"^ 

<a/3> 



(u^ - m1) • t"^ 
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<a/37>ni 

^<5^fe^«-«l). (A 9) 

■^"^<a/37>v 

The supplementary factor 1/V3 which does not appear in (2.15) is due to the 
factor \/3 between the equilibrium distances P — a and /? — 7 (where /3 is the apex 
of the angle). The coefficients of the expansion are e"'' = fjg'^ + ^6^^ - i^g'^, 
e"-P = -^6f - 5f + ISf, and k„/e = ^'(-662 + 1266). 

(ii) Invariance requirements 

An infinitesimal rotation around an axis parallel to n of the undistorted lattice 
yields u" = 9n x r", parallel to V. Hence, it" — it'' — —6£n x for nearest 
neighbours, and u^ — u^ = —O^V^nxf^f^ for next-nearest neighbours. Substituting 
these expressions into (2.20) yields F^^)/{e'^e) = -|V3iV2(G^ + 0% + 3Gg) + 
3iV(4G + \/3G| + 3\/3Gg) - 3(4G + \/3G^). Setting F^^^ = for all N gives 
V3G^ = 8G, V^G^ = 4G and V3Gg = -4G. Finally, imposing the same rotation 
on a properly chosen initial state allows to deduce that G must be set equal to zero 
to satisfy the rotational invariance. 
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